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Abstract 



| Correlation of charm quark - charm antiquark in 7p scattering 

are calculated in the fcj-factorization approach. We apply different 
unintegrated gluon distributions (uGDF) used in the literature. The 
\ results of our calculations are compared with very recent experimental 

results from the FOCUS collaboration. The CCFM uGDF developed 
recently by Kwiecihski et al. gives a good description of the data. 
New observables are suggested for future studies. Predictions and 
perspectives for the HERA energies are presented. 

In recent years a lot of activity was devoted to the description of the 
photon-proton total cross section (or F 2 structure function) in terms of the 
unintegrated gluon distribution functions (uGDF) (see e.g. [U 12] and refer- 
ences therein). In some of the analyses also inclusive charm quark (or me- 
son) were considered [3|. Although the formalism of uGDF is well suited for 
stydying more exclusive observables, only very few selected cases were con- 
sidered in the literature. A special example is azimuthal jet-jet correlations 
in photon-proton scattering |U 03 E] . In this case one samples in a nontrivial 
way simultaneously the x and k t dependences of uGDF. The HI collaboration 
at HERA has measured very recently such correlations [Zj. Very similar anal- 
ysis was performed in the past also for open charm production at somewhat 
lower energies jH]. Recently the FOCUS collaboration at Fermilab provided 
new precise data for charm-anticharm correlations jHj- It is our aim here to 
analyze the charm-anticharm correlations in terms of uGDF. In the present 
paper we wish to compare results for different uGDF available in the litera- 
ture. While the total cross section depends on small values of x, the high-p f 
jets and/or heavy quark production test gluon distributions at somewhat 
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larger x. Only some approches from the literature are applicable in this re- 
gion. In particular, we wish to test results based on CCFM unintegrated 
parton distributions developed recently by Kwiecihski et al.[10[ ITT]. 

The total cross section for quark-antiquark production in the reaction 
7+p— > Q + Q + X can be written as [Sj 
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In the formula above f g (x, k 2 ) is the unintegrated gluon distribution with the 
convention from Ref. 6j. 1 The gluon transverse momentum is related to the 
quark/antiquark transverse momenta pij and p2,t as: 
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The unintegrated gluon distribution f g is evaluted at 
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It is obvious that at larger transverse momenta of quarks and/or heavy quark- 
antiquark production one samples larger values of x g than in the case of the 
total photon-proton cross section. 

The choice of (j% is not essential for the discussion in the present paper and 
will be discussed elsewhere. In the present calculation the scale of running 
coupling constant in @ is taken to be k 2 and the freezing prescription from 
[T2*] is used. The latter prescriptions are not very important when study- 
ing correlations. They may be important, however, for the integrated cross 
sections, to be studied elsewhere [T3] . 

Thus the basic ingredient of our approach are unintegrated gluon dis- 
tributions. Different models of uGDF have been proposed in the literature 

1 Two different conventions are used throughout the literature 
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(see for instance Pfl |2j and references therein). The main effort has been 
concentrated on the small-x region. While the total cross section is the gen- 
uine small-x phenomenon (x < 10 -3 ), the production of charm and bottom 
quarks samples rather the intermediate-x region (x~ 10~ 2 - 10 _1 ) even at 
the largest available energies at HERA. It is not obvious a priori if the meth- 
ods used are appropriate for the intermediate values of x. In the present 
approach we shall present results for a few selected gluon distributions from 
the literature. For illustration we shall consider the simple BFKL [03, the 
saturation model used to study HERA photon-proton total cross sections 
[16J (GBW), and the saturation model being often used recently to calculate 
particle production in hadron-hadron collisions (KL). These three model 
approaches are expected to be valid for small, not very well specified, values 
of x. At somewhat larger values of x all these models are expected to break. 
This may happen already at the FOCUS energy W = 18.4 GeV. Clearly 
an extrapolation may be needed. As in Ref . [Tlj one can try to extend the 
applicability of these small-x models by multiplying the model distributions 
by a phenomenological factor (1 — x) n . In principle, the value of n could be 
adjusted to inclusive spectra at lower energies. The choice of n is, however, 
marginal for the correlations studied in the present paper. We shall not dis- 
cuss here the details of the different approaches. A more detailed discussion 
can be found in Ref. ^1] . 

In addition, we shall consider two other approaches adequate for intermediate- 
x region. The CCFM approach seems to be the best tailored for this purpose. 
It was shown in Refs. [TO} ITT] how to solve the one-loop CCFM equation in 
the impact parameter representation. 

The unintegrated parton distributions used in the present paper were ob- 
tained by solving the Kwiecihski CCFM equations jTUJ HJj using LO GRV98 
collinear distributions 22 j as the input for the evolution. By construction 
this procedure assures that our uPDF provide a good description of the Fi 
structure function data. 

The solution of the CCFM equation depends on three variables f g = 
f g (x,b, fi 2 ). 2 The familiar momentum representation unintegrated gluon 
distribution can be obtained via Fourier-Bessel transform 

f q (x,K 2 ,fi 2 ) = J exp (iKb) f q (x,b,fi 2 ) d 2 b . (6) 

As already mentioned in the introduction, it is our intention here to use 
uGDFs fg CFM {x, b, fj 2 ) which fulfill the b-space (one- loop) CCFM equations 
jTUl E]. However, the perturbative solution fg CFM (x, b, /i 2 ) does not include 

2 In the present paper we shall use the notation / instead of / as in Refs. |18M19| . 
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nonperturbative effects such as, for instance, intrinsic momentum distribu- 
tion of partons in colliding hadrons. In order to include such effects we 
propose to modify the perturbative solution f GCFM (x,b, fj 2 ) and write the 
modified gluon distribution f g (x,b,fi 2 ) in the simple factorized form 

fg( x ,b, //) = fg CFM (x, b, //) • F^ p (b) . (7) 

In Ref. |2Uj. two different functional forms for the nonperturbative form factor 

F NP {h) = F NP {h) = exp or exp ^_£) (g) 

identical for all species of partons were used. In Eq.(JHJ) b (or b e ) is the only 
free parameter. The parameters were roughly adjusted in [20 to describe 
production of W and Z bosons in nucleon-nucleon collisions. In the present 
note we shall show only results for the Gaussian form factor. The dependence 
on the choice of the form factor will be studied elsewhere. 

The resummation formulae [21] and the unintegrated parton distribution 
formulae for Higgs [19J and gauge boson [2U] have identical structure if the 
following formal assignment is made: 

~f g GR {x,b^ 2 ) = ^F g NP (fi,b,x)[xg(x 1 ^(b)) + ..]exp(^S g (b, f i)\ , 
f^ GR (x,b,^) = ^F q NP (^b,x)[xq(x 1 ^(b)) + ...}exp(^S q (b. if ,) 
f^ R (x,b,v 2 ) = h q r^,b,x)[xq(x, f ,(b)) + ..]exp(^S, / (b,^ . (9) 

The index SGR above stands for "soft-gluon resummation". The explicit 
expressions for S g , S q and Sq can be found in Ref. [TH} EHj . where in addi- 
tion similarities and differences between Kwiecihski CCFM and soft gluon 
resummation are discussed. 

The ^-dependent unintegrated distributions of gluons corresponding to 
the b-space resummation can be then obtained through the Fourier-Bessel 
transform 

f^ GR (x, K 2 ,Q 2 ) = J dbb J (K,b)fg GR (x, b, Q 2 ). (10) 

With the simple Ansatz (JSJ) for F^ 9 the whole scale [i 2 dependence resides 
exclusively in the Sudakov-like form factor. For brevity, we shall call the 
gluon distribution in Eq. ()10|) the "resummation gluon distribution" . 

Before we go to charm quark - charm antiquark correlations we wish 
to show the results for inclusive spectra of c or c. In Fig^ we show the 
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distributions for different uGDF for low (W=18.4 GeV, left panel) and 

high (W=200 GeV, right panel) energies. At the lower energy the GBW 
uGDF gives somewhat steeper pf distribution than the other uGDF. At the 
high energy the slope of the p\ distributions decreases. Quite similar slopes 
are obtained for different uGDF. Therefore this observable is not the best one 
to test models of uGDF. We shall show below that more exclusive correlation 
observables are more sensitive tests of models/parametrizations of uGDF. 

In the present analysis we do not put any restrictions on heavy quark or 
heavy antiquark transverse momenta pij and p2,t- A detailed analysis of the 
effect of such cuts on the results will be presented elsewhere [T5j . 

The azimuthal correlation functions w(<j>) defined as: 



S%W) 
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da 



w(d>;W)= r Z T J J2 , (11) 



W) = J dpi J dpi J dz . a(W,p l!t ,p 2>t ,z) (12) 

in order that 

r w{<p]W) d(f> = 1 (13) 



for two energies ofW = 18.4 GeV (FOCUS) and W = 200 GeV (HERA) are 
shown in Fig0 The GBW-glue (thin dashed) gives too strong back-to-back 
correlations for the lower energy. Another saturation model (KL, [TJj) pro- 
vides more angular decorrelation, in better agreement with the experimental 
data. The BFKL-glue (dash-dotted) provides very good description of the 
data. The same is true for the CCFM-glue (thick solid) and resummation- 
glue (thin solid). The latter two models are more adequate for the lower 
energy. The renormalized azimuthal correlation function ()13|) for BFKL, 
GBW and KL models are almost independent of the power n in extrapo- 
lating to larger values of x g . In the present paper, in calculating the cross 
section with the Kwiecihski CCFM uGDF for simplicity we have fixed the 
scale for fx 2 = 4 m 2 . Allowing for dependence of the scale fi 2 on kinematical 
variables such as x or k 2 would make the calculation very time consuming. 
The sensitivity to the choice of the scale will be discussed in detail elsewhere. 
For comparison in panel (b) we present predictions for W = 200 GeV. Ex- 
cept of the GBW model, there is only a small increase of decorrelation when 
going from the lower fixed-target energy region to the higher collider-energy 
region. 
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In the present paper we completely ignore the resolved photon component 
j2H|. The latter should be, however, negligible for the FOCUS fixed target 
experiment [Oj. 

Not only azimuthal correlations are interesting. In general, the integrand 
of Eq.((H) depends on four independent kinematical variables 4>,p 2 t ,pl, z 
(other combinations of the kinematical variables are also possible). In par- 
ticular, the formula ((TJ can be rewritten in the form 

a^QQ(W) = J dpi J dpi w(pl,pl- W) (14) 

where the two-dimensional correlation function 

w(pl,pl;W) = J d(f> dz fa{x ^ K2) .a(W,p ht ,p 2it ,z). (15) 

In FigOJwe present some examples of w{p\ t , p\ t ) at W=18.4 GeV for different 
uGDF. The maps for different uGDF differ in details. The distribution for 
the GBW gluon distribution is concentrated along the diagonal p\ t = pi 
and in this respect resembles the familiar collinear leading-order result. 3 
The other three distributions have a sizeable strength at the phase-space 
borders for p\ t f» or p\ t « 0. The KL gluon distribution gives in addition 
some enhancement at p\ t ~ pi, especially at large transverse momenta. 
Experimental studies of such maps could open an interesting new possibility 
to test models of uGDF in a more detailed differential fashion. In principle, 
such studies will be possible with HERA II runs at DESY. The leading-order 
approach of Kwiecihski contains the higher-order corrections via evolution 
equations. However, in contrast to NLO collinear approach it can be applied 
even in the region p\^ t = P2,t- 

At present experimental luminosities (statistics) one may have a problem 
to explore the whole two-dimensional maps shown in FigJHJ In this case a 
more global variable could be useful. In order to quantify the spread over 
Pit x P%t plane and/or departure from the diagonal (LO collinear approach) 
we propose a new variable which can be interpreted as a measure of deviations 
from the equal-length momenta defined as: 

f(pL, > k P 2 min , W) = a{P ™ >^ W) , (16) 

where p 2 max = max(pf >t ,p| t ) and p 2 min = min(p^,p^). For example with k 
= 2 one obtains respectively: GBW: 0.01, KL: 0.45, BFKL: 0.63, K(CCFM): 

3 It is well known that the collinear NLO calculation is not reliable for pi t = p2,t- 
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0.60, resum: 0.33. The quantity f{j? max > W) is sensitive to both 

perturbative 7^ p 2 ,t) and nonperturbative {p^ t ~ P2,t) processes and 
therefore reflects their interplay. We believe that the FOCUS collaboration 
could reprocess their present data in order to obtain analogous fractions for 
their {p 2 £>iP 2 jj) distributions. 

In the leading-order collinear approach (without parton showers included) 
the transverse momenta of two jets add up to zero. It is not the case for our 
leading-order ^-factorization approach. In FigE] we present normalized to 
unity distribution in p 2 +) where p + = p\ + p2- Due to momentum conser- 
vation, in our approach the sum of transverse momenta is directly equal to 
transverse momentum of the gluon (p + = R). This means that the distribu- 
tion in p 2 + directly probes the transverse momentum distribution of gluons. 
The situation here is very similar to the situation for the azimuthal angle 
distribution. The GBW gluon distribution which discribes very well the to- 
tal 7* — p cross section gives very steep distribution in p 2 + in comparison to 
other uGDF. We expect that the inclusion of QCD evolution effects, like in 
Ref . [21] for instance, should change this result and lead to somewhat broader 
distributions. The Kwiecihski CCFM gluon distribution gives the best de- 
scription of the FOCUS data jHj. We expect that this approach is suitable 
for x > 0.01. Although the other models give also reasonable description of 
the FOCUS w(p 2 + ) data, one should remember that their application for the 
low-energy data is somewhat unsure. 

For completeness in FigJ5]we present normalized to unity distributions in 
the square of p_ = p\ — p2- Here, however, the differences between different 
uGDF are smaller than in the previous two cases. Nevertheless the corre- 
sponding data would be a new possible observable to verify the unintegrated 
gluon distributions. 

In summary, we have shown that the analysis of kinematical correlations 
of charm quarks and antiquarks opens new possibilities for verifying models of 
uGDF. The recently measured data of the FOCUS collaboration at Fermilab 
allows one to study the unintegrated gluon distribution in the intermediate-x 
region. Many models of uGDF used in the literature are constructed rather 
for small values of x and its application in the region of somewhat larger x (x 
> 0.05) is questionable. The unintegrated gluon (parton) distribution which 
fulfill the CCFM equations, developed recently by Kwiecihski et al., describe 
the data fairly well. It can be expected that the correlation data from the 
HERA II runs will give a new possibility to verify the different models of 
unintegrated gluon distributions in a more detailed way. 
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Figure 1: Inclusive distribution of charm quarks /antiquarks as a function 
of the corresponding transverse momentum. We present results for GBW 
(thin dashed), KL (thick dashed), BFKL (dash-dotted) and Kwiecinski CCFM 
(solid). 




Figure 2: Azimuthal correlations between c and c. The theoretical results 
are compared to the recent results from ^ (fully reconstructed pairs). The 
notation is the same as in FigU\ 
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Figure 3: Some examples of the two-dimensional maps w(p\ t ,p\ t ) for W = 
18.4 GeV: a) GBW, b) KL, c) BFKL, d) CCFM, e) resummation. 
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Figure 4: p 2 + distribution ofc — c. The theoretical results are compared to the 
recent results from [9] (fully reconstruced pairs). The notation is the same 
as in FigU\ 
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